Strong Dependence of Ultracold Chemical Rates on Electric Dipole Moments 
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We use the quantum threshold laws combined with a classical capture model to provide an analyt- 
ical estimate of the chemical quenching cross sections and rate coefficients of two colliding particles 
at ultralow temperatures. We apply this quantum threshold model (QT model) to indistinguishable 
fermionic polar molecules in an electric field. At ultracold temperatures and in weak electric fields, 
the cross sections and rate coefficients depend only weakly on the electric dipole moment d induced 
by the electric field. In stronger electric fields, the quenching processes scale as d 4t ^ L+ ^ where 
L > is the orbital angular momentum quantum number between the two colliding particles. For 
p— wave collisions (L = 1) of indistinguishable fermionic polar molecules at ultracold temperatures, 
the quenching rate thus scales as d 6 . We also apply this model to pure two dimensional collisions 
and find that chemical rates vanish as d~ 4 for ultracold indistinguishable fermions. This model 
provides a quick and intuitive way to estimate chemical rate coefficients of reactions occuring with 
high probability. 



I. INTRODUCTION 

Ultracold samples of bi-alkali polar molecules have 
been created very recently in their ground electronic 1 E, 
vibrational v = 0, and rotational N = states 
This is a promising step before achieving Bose-Einstein 
condensates or degenerate Fermi gases of polar molecules, 
provided that further evaporative cooling is efficient. For 
this purpose, elastic collision rates must be much faster 
than inelastic quenching rates. This issue is somewhat 
problematic for the bi-alkali molecules recently created, 
since they are subject to quenching via chemical reac- 
tions. If a reaction should occur, the products are no 
longer trapped. 

For alkali dimcrs that possess electric dipole moments, 
elastic scattering appears to be quite favorable, since elas- 
tic scattering rates are expected to scale with the fourth 
power of the dipole moment 0, [H ■ Inelastic collisions of 
polar species can originate from two distinct sources. The 
long-range dipole-dipole interaction itself is anisotropic 
and can cause dipole orientations to be lost. This kind 
of loss generally leads to high inelastic rates, and is re- 
garded as the reason why electrostatic trapping of polar 
molecules is likely not feasible @. Moreover, these colli- 
sions also scale as the fourth power of dipole moment in 
the ultracold limit [|| , meaning that the ratio of elastic to 
inelastic rates does not in general improve at higher elec- 
tric fields. This sort of loss can be prevented by trapping 
the molecules in optical dipole traps. 

More serious is the possibility that collisions are 
quenched by chemical reactions. Chemical reaction rates 
are known to be potentially quite high even at ultra- 
cold temperatures ■ Indeed, for collision ener- 
gies above the Bethe-Wigner threshold regime, it ap- 
pears that many quenching rates, chemical or other- 
wise, of barrier less systems are well described by apply- 
ing Langevin's classical model [l8|]. In this model the 
molecules must surmount a centrifugal barrier to pass 
close enough to react, but are assumed to react with unit 
probability when they do so. This model has adequately 



described several cold molecule quantum dynamics cal- 
culations 0, EH, E El [13 ■ 

Within the Bethe-Wigner limit, scattering can be de- 
scribed by an elegant Quantum Defect Theory (QDT) 
approach [l^-[22j ■ This approach makes explicit the dom- 
inant role of long-range forces in controlling how likely 
the molecules are to approach close to one another. Con- 
sequently, quenching rate constants can often be written 
in an analytic form that contains a small number of pa- 
rameters that characterize short-range physics such as 
chemical reaction probability. For processes in which the 
quenching probability is close to unity, the QDT th eory 
provides remarkably accurate quenching rates p3l . |24| . 
For dipoles, however, the full QDT theory remains to be 
formulated. 

In this article we combine two powerful ideas - sup- 
pression of collisions due to long-range physics, and high- 
probability quenching inelastic collisions for those that 
are not suppressed - to derive simple estimates for in- 
elastic/reactive scattering rates for ultracold fermionic 
dipoles. The theory arrives at remarkably simple ex- 
pressions of collision rates, without the need for the full 
machinery of close-coupling calculations. Strikingly, the 
model shows that quenching collisions scale as the sixth 
power of the dipole moment for ultracold p— wave colli- 
sions. On the one hand, this implies a tremendous de- 
gree of control over chemical reactions by simply varying 
an electric field, complementing alternative proposals for 



electric field control of molecule-molecule [25| or atom- 
molecule [26[ chemistry. On the other hand, it also im- 
plies that evaporative cooling of polar molecules may be- 
come more difficult as the field is increased. In section II, 
we formulate the theoretical model for three dimensional 
collisions. In section III, we apply this model to pure 
two dimensional collisions and conclude in section IV. In 
the following, quantities are expressed in S.I. units, un- 
less explicitly stated otherwise. Atomic units (a.u.) are 
obtained by setting h = Atteq = 1. 
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II. COLLISIONS IN THREE DIMENSIONS 



A. Cross sections and collision rates 

In quantum mechanics, the quenching cross section of 
a pair of colliding molecules (or any particles) of reduced 
mass /i for a given collision energy E c and a partial wave 
L, Ml is given by 



a 



qu 

L,Ml 



2\iE c 



I?TmJ 2 xa 



(1) 



where T qu is the transition matrix element of the quench- 
ing process, \TI u Ml \ 2 represents the quenching proba- 
bility, and the factor A represent s sy mmetrization re- 
quirements for identical particles [27| . If the two col- 
liding molecules are in different internal quantum states 
(distinguishable molecules), A = 1 and if the two 
colliding molecules are in the same internal quantum 
state (indistinguishable molecules), A = 2. The to- 
tal quenching cross section of a pair of molecules is 
a qu _ & q LM L ■ The quenching rate coefficient of a 

pair of molecules for a given temperature T (collisional 
event rate) is given by 
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where 



f{v) =4tt 



2TTk B T 
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v 2 exp[-{fiv 2 )/{2k B T)} (3) 



is the Maxwell-Boltzmann distribution for the relative 
velocities for a given temperature and k B is the Maxwell- 
Boltzmann constant. The total quenching rate coeffi- 
cient of a pair of molecules is K qu = J2 L M K'£ 1 Ml . 
To avoid confusion, we will also write the correspond- 
ing rate equation for collisions between distinguishable 
and indistinguishable molecules. First, we consider colli- 
sions between two distinguishable molecules in quantum 
states a and b (A = 1 in Eq. ([I])). During a time dt = t, 
where r is the time of a quenching collisional event, the 
number of molecules N a lost in each collision is one and 
the number of molecules N b lost in each collision is one. 
Then dN a /dt = -1/r and dN b /dt = -1/r. The volume 
per colliding pairs of molecules is V/(N a N b ), where V 
stands for the volume of the gas. During the time r, the 
quenching collisional event is associated with a volume 
< a qu xv > xt = K qu xr. By definition of t, this volume 
should be equal to the one occupied by just one colliding 



pair of molecules. Then we get K 



qu 



x r 



V/{N a N b ). 



The rate equation for the number of molecule N a or N b 
is then given by 



dN, 



= -K 



NgN b 

V ' 



< 4 » 

If n a = N a /V and n b = N b /V arc the densities of 
molecule a and b in the gas, then 

dn„ 



dt 



_ K qu 



x n a n b . 



(5) 




FIG. 1: (Color online) Effective potential barrier V(R) as a 
function of the intermolecular separation R. Vb and R b denote 
the height and the position of the centrifugal barrier. 



We consider now the case of collisions between two in- 
distinguishable molecules (A = 2 in Eq. ([1])). During 
the time dt = t the number of molecules N lost in each 
collision is two. Then we get dN/dt = — 2 jr. The vol- 
ume per colliding pairs of molecules is V/{N(N — I)/2) 
where we have taken into account the indistinguisha- 
bility of the molecules. For the same reason explained 
above, the volume associated with the collisional event 
during the time r should be equal to the volume occu- 
pied by just one colliding pair of molecules. And then 
K qu xr = V/(N(N - l)/2). The rate equation for the 
number of molecule N is then given by 



dN _ N(N -l)/2 

dt V 



If n = N/V and N(N - I) w TV 2 , then 



^ = -K qu x n 2 
dt 



(6) 



(7) 



B. Quantum threshold model 

We consider the case of two identical ultracold 
fermionic polar molecules, as has been achieved very re- 
cently for KRb dimers @, ,28] in their ro-vibronic v — 
0, N — 0) ground state. Under these circumstances, be- 
cause of Fermi exchange symmetry, the relative orbital 
angular momentum quantum number L between the two 
molecules must take odd values L = 1,3,5,7.... These 
molecules are polar molecules and can be controlled by 
an electric field £. In the usual basis set of partial waves 
\LMl), the long-range behavior of two colliding polar 
molecules in a presence of an electric field is governed by 
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an interaction potential matrix whose elements are 

(LM L \V(R)\L'M' L ) = 

h 2 L(L + l) C 6 \ x 

C 3 (L,L';M L ) 



R 3 



■>M L ,M'r 



(8) 



where R denotes the distance between the two molecules. 
The diagonal elements represent effective potentials 
for the colliding molecules and the non-diagonal el- 
ements represent couplings between them. The co- 
efficient Cq is the van der Waals coefficient, as- 
sumed to be isotropic in the present treatment. The 
—C3/R 3 is the term corresponding to the electric dipole- 
dipole interaction expressed in the partial wave ba- 
sis set (LML\Vdd(R,Q,<p)\L'M' L ) between two polar- 
ized molecules in the electric field direction, with 
V dd (R, 9, ip) = d 2 (l - 3 cos 2 6<)/(47r£o R 3 ), where d = d(£) 
is the induced electric dipole moment, and 9, tp represent 
the relative orientation between the molecules. In the 
basis set of partial waves, C3 takes the form 

j2 

C 3 (L,L';M L ) = a(L,L';M L 



47T£ 

2 (-l) Mt V2L + 1 V2L' + 1 
L 2 L' 
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The large bracket symbols denote the usual 3— j coeffi- 
cients. The coefficient a is introduced to simplify further 
notations. The combination between repulsive and at- 
tractive terms in the effective potentials (diagonal terms) 
of Eq. ((8]) generate a potential barrier of height Vb which 
is plotted schematically in Fig. [TJ The height of this bar- 
rier plays a crucial role as it can prevent the molecule 
from accessing the short range region where reactive 
chemistry occurs. 

The quantum threshold (QT) model consists of two 
conditions. First, for E c < Vb, we use the Bethe-Wigner 
threshold laws [29|, Ho| for ultracold scattering. Second, 
we use the classical capture model (Langevin model) [l8| 
to estimate the probability of quenching for E c > Vb. A 
classical capture model is indiferent to collision energies 
E c < Vb since the barrier prevents the molecules from 
coming close together. In real-life quantum scattering, 
collisions do occur at these energies due to quantum tun- 
neling, and they are the ones relevant to ultracold colli- 
sions. Moreover, collisions in this energy regime are dic- 
tated by the the Bethe-Wigner quantum threshold laws. 
For quenching collisions, the threshold laws [29M3l1 | state 

that iTjfjtf I oc E c 2 . For E c > Vt, a classical capture 
model will guarantee to deliver the molecule pair to small 
values of R, where chemical reactions are likely to occur 
with unit probability (see Fig. [I]). Following this classi- 
cal argument, we will assume that when E c > Vb, the 
quenching probability reaches unitarity |T 9 "| 2 = 1. Us- 
ing this assumption together with the quantum threshold 



laws, the QT quenching tunneling probability below the 
barrier can be written as 



L,M L \ 



L + l/2 



(10) 



Consequently, the quenching cross section and rate coef- 
ficient are approximated by 



H 2 TT 



'LMl 



K qu 
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(11) 



for E c < Vb- The QT model has the simple and intuitive 
advantage of showing how the cross sections and rate co- 
efficients scale with the height of the entrance centrifugal 
barrier. For E c > Vb, it is easy to find the correspond- 
ing expression of the cross section in Eq. ([T]) by setting 

| T g«|2 = L The crosg section 

®1lMl w ^ reacrl the uni- 
tarity limit at E c > Vb- It is also easy to find the cor- 
responding expression of the rate coefficient in Eq. (|2j). 
The QT model is general for any collision between two 
particles provided that there is a barrier in the entrance 
collision channel and that chemical reactions occur with 
near unit probability at short range. The only informa- 
tion on short range chemistry is that chemical reactions 
occur at full and unit probability and the only informa- 
tion on long range physics is provided by the height of the 
entrance barrier Vb- The QT model describes the back- 
ground scattering process, it does not take into account 
scattering resonances. Note that the model will not be 
appropriate in the present form for barricrlcss (s— wave) 
collisions since Vb = 0. For this particular type of colli- 
sions that do not possess a centrifugal barrier, the QDT 
theory can be usefully applied [U [H, HH . The present 
form of the QT model does not take into account the 
anisotropy of the intermolecular potential at intermedi- 
ate range and/or the electronic and nuclear spin structure 
of the molecular complex but remains suitable as far as 
the entrance centrifugal barrier takes place at long range. 
The QT model will have to be modified if longer range 
interactions takes place. For example, collisions between 
N = and N = 1 polar molecules can have long range 
interactions between hyperfine states due to dipolar and 
hyperfine couplings [23, [3^. However, for collisions be- 
tween rotationless TV = polar molecules, the hyperfine 
couplings are weak and the QT model can be applied 
without further modifications. 



C. Rates in zero electric field 

In the absence of an electric field in Eq. (JSJ) , the long 
range potential reduces to a diagonal term in the basis set 
of partial waves. The position and height of the barrier 
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FIG. 2: (Color online) Quenching cross section of 39 K + 
39 K2 as a function of the collision energy for the partial wave 
L = 1 — 5: (i) calculated with a full quantum calculation (solid 
lines), reproduced from Ref. [l0|] (ii) using the QT model 
(dashed lines) (iii) fitting the QT model (dotted line), using 
p = 0.34 in Eq. (fH)) . 
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FIG. 3: (Color online) Quenching probability of 39 K + 39 K 2 
as a function of the collision energy for the partial wave L = 
I: comparison between the full quantum calculation (solid 
line) and the QT model (dashed line). The fitted QT model 
appears as a dotted line. The height of the barrier Vb and the 
corrected height 7 X Vf, (7 = 2.06) appear as vertical lines. 



are given by 



Rb — 




(12) 



We can insert Eq. (|12[) in Eq. (|TTj) to get analytical forms 
of the quenching cross section or rate coefficient. For two 
indistinguishable fcrmionic polar molecules at ultracold 



temperatures when L = 1, we get 



K qu 



it ^3 1 _VC|N 1/4 



k B T x A. (13) 



In Eq. (|TT]), we used the fact that < E c >= 3k B T/2 in 
three dimensions. Note that to get the overall contribu- 
tion for a given L, we have to multiply Eq. (|13[) by the 
degeneracy factor (2L + 1) corresponding to all values of 
Ml ■ We can get similar expressions for any partial wave 
L. 

To test the validity of the model, we compare in Fig. [2] 
the quenching cross sections of 39 K( 2 5*) + 39 Fv2( 3 £+, v = 
1,N = 0) collisions as a function of the collision en- 
ergy for the partial waves L = 1 — 5: (i) calculated in 
Ref. [lfjj with a full quantum time-independent close- 
coupling calculation based on hypcrspherical democratic 
coordinates (HJ and the full potential energy surface of 
K3 (solid lines) (ii) using the simple QT model (dashed 
lines) with a value of C§ = 9050 a.u. given in Ref. [l(| 
(1 a.u. = 1 Eh cto where Eh is the Hartrcc energy and 
an is the Bohr radius). In this example, the QT model 
provides an upper limit to the cross sections. This is 
due to the fact that the quenching cross section does not 
reach a maximum value at the height of the barrier Vb, 
but rather at somewhat higher energy, say 7 x Vb, with 
7 > 1 (see Ref. [l(|). For all partial waves, there is a 
worse agreement for collision energies in the vicinity of 
the height of the barrier where the passage from the ul- 
tralow regime to the unitarity limit is smoother than for 
the QT model. This smoother passage is visible in Fig. m 
for the full quantum L = 1 quenching probability com- 
pared to the QT model, which has a sharp corner in the 
vicinity of Vb- 

To account for more flexibility in the QT model, one 
can replace Vj, in Eq. (fTUf by 7 x Vb (7 > 1), and use the 
coefficient 7 as a fitting parameter to reproduce either 
full quantum calculations or experimental observed data. 
Alternatively, we can correct the QT quenching tunneling 
probability with an overall factor p, 



\rpqu 
\ ± L.M 1 



I 2 

lit 



[ v h 



L+l/2 



(14) 



with p = ry-(- L+L / 2 ) . p < 1 can be interpreted as the 
quenching probability reached at the height of the bar- 
rier Vb in the QT model, rather than the rough full 
unit probability (p = 1). As an example, we find that 
7 w 2.06 reproduces the quantum L = 1 partial wave 
cross section for 39 K + 39 K2 (dotted line in Fig. [2] and 
Fig.[3|). This yields a maximum quenching probability of 



P = 7 



-3/2 



0.34 instead of 1. In other words, the QT 



i 2.96 higher than the full 
39 K2 collisions at ultralow 



model is only a factor of p~ x 
quantum calculation for 39 K 
energies. 

Given the fact that full quantum calculations are com- 
putationally demanding @-[Tj] and impossible at the 
present time for alkali molecule-molecule collisions, the 
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accuracy of the QT model is satisfactory and can be a 
quick and powerful alternative way to estimate orders of 
magnitude for the scattering observables. Besides, agree- 
ment between the QT model with experimental data or 
full quantum calculations is expected to be satisfactory 
for collisions involving alkali species, because it is likely 
that short range quenching couplings will dominate and 
lead to high quenching probability in the region where 
the two particles are close together [l6j]. Very recently, 
Eq. (fT3"| of the QT model has been applied for the eval- 
uation of ultracold chemical quenching rate of collisions 
of two 40 K 87 Rb molecules in the same internal quantum 
state, and provided good agreement with the experimen- 
tal data 1321. 
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D. Rates in non-zero electric field 

1. Numerical expressions 



KRb( 1 £,v=0,N=0) + KRb( 1 E,v=0,N=0) ; L=1 




d(D) 

FIG. 4: (Color online) Diabatic (dashed lines) and adiabatic 
(solid lines) barrier heights V b ' a as a function of the induced 
dipole moment d for the partial waves L = 1,Ml = (red 
curves) and L = 1, \Ml\ = 1 (blue curves). 

In the presence of an electric field in Eq. (|8|) , the long- 
range interaction potential matrix is no more diagonal 
and couplings between different values of L occur. Ml is 
still a good quantum number. A first approximation (di- 
abatic approximation) consists of neglecting these cou- 
plings and using only the diagonal elements of the di- 
abatic matrix directly. Then one can find numerically 
for which R the centrifugal barriers are maximum and 
evaluate the height of the diabatic barriers Vjf. This is 
repeated for all values of the induced dipole moment d. 
A second approximation (adiabatic approximation) is to 
diagonalize this matrix (including the non-diagonal cou- 
pling terms) for each R and again find the maximum of 
the centrifugal barriers to get the height of the adiabatic 
barriers V b a . As an example, we compute these barrier 
heights for 40 K 87 Rb— 40 K 87 Rb collisions, using a value 



FIG. 5: (Color online) Quenching rate coefficients of two in- 
distinguishable fermionic polar 40 K 87 Rb molecules as a func- 
tion of the induced electric dipole moment for L = 1 and 
for a temperature of T = 350 nK (black curves). The rates 
have been calculated using the barrier heights of Fig. [4] The 
red lines represent the L = 1,Ml ~ partial wave contri- 
bution. The blue lines represent the sum of L = 1,Ml = 1 
and L = 1, Ml = — 1 partial wave contributions. The dashed 
lines represent the rates calculated with the diabatic barriers 
while the solid lines with the adiabatic barriers (see text for 
detail). The total, Ml ~ and \Ml\ = 1 curves have been 
indicated in the left hand side. 



of C 6 = 16130 a.u. [H. We P lot in Fig. g] the heights 
of the diabatic (dashed lines) and adiabatic (solid lines) 
barriers for the quantum numbers Ml = (red curves) 
and | Ml | = I (blue curves). The adiabatic barriers have 
been calculated using five partial waves L = 1 — 9 in 
Eq. <JSj> - The effect of the couplings can be clearly seen in 
this figure by comparing diabatic and adiabatic barriers. 
Especially for the \M L \ = 1 case for d « 0.16 D (1 D 
= 1 Debyc = 3.336 10 -30 Cm), couplings with higher 
partial waves make the adiabatic barrier decrease as the 
dipole increases while the diabatic barrier continues to 
increase. 

Using these heights of the barriers, we use Eq. (jTTJ) 
to plot in Fig. [5] the total quenching rate coefficients 
(black curves) as a function of d for two indistinguishable 
fermionic polar 40 K 87 Rb molecules in the same quantum 
state for L = 1 and at a typical experimental temperature 
of T = 350 nK [H. For T = 350 nK, the mean collision 
energy < E c >= 3ksT/2 = 525 nK, and the maximum 
dipole moment for which Vb < 525 nK (that is for which 
Eq. (jTTJ) does not apply anymore) is around d « 0.24 D 
(see Fig. [4j. The dashed curves correspond to rates cal- 
culated with the diabatic approximation while the solid 
curves correspond to rates calculated with the adiabatic 
approximation. The Ml = contribution is plotted in 
red and the contribution of Ml = +1 and Ml = — I is 
plotted in blue. The rates highly reflect the behavior of 
the centrifugal barriers in the entrance collision channel. 
When the barrier increases with the dipole, it prevents 
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FIG. 6: (Color online) Same as Fig. [5] but we use the an- 
alytical expressions for the rates (see text for detail). The 
total, Ml = and \Ml\ = 1 curves have been indicated in 
the left hand side. The individual analytical curves have been 
indicated in the right hand side by roman numbers. 



the molecules from getting close together and the quench- 
ing rates decreases. When the barrier decreases, the tun- 
neling probability is increased allowing the molecules to 
get close together, and the quenching rates increases. 



2. Analytical expressions 

In order to have an intuitive sense of how the chemi- 
cal quenching rate scales with the induced dipole moment 
(and the electric field) , we evaluate analytical expressions 
of the barriers and the rates as it has been done in the 
previous section for a zero electric field. The analyti- 
cal expression of the height of the diabatic barrier V b is 
complicated by the occurrence of two distinct long-range 
potentials in the diagonal matrix term of Eq. ([8]). We cir- 
cumvent this difficulty by looking in the two limits where 
one dominates over the other. For small electric fields, we 
use the zero electric field limit discussed in the preceding 
section by setting C3 = 0. For larger electric fields we 
ignore the Cq coefficient in Eq. ((8|) if the electric dipole- 
dipole interaction is attractive (positive C3). We ignore 
the centrifugal term in Eq. ((5J) if the electric dipole-dipolc 
interaction is repulsive (negative C3). These two cases 
are discussed below. In between, to accommodate the 
transition between the low-field and high-field limit, we 
will simply add the rate coefficients derived in the two 
limiting cases. 

For positive C3 coefficients in Eq. ([9]), —C3/R 3 is at- 
tractive in Eq. ([8]). For L = 1 partial waves for example, 
this occurs when Ml = 0, which favors an attractive 
orientation of dipoles. We consider 



R 3 



> 



R 6 



(15) 



in Eq. (JH). In this case, the position and height of the 
barrier are given by 



Rb 



3/iC 3 



h 2 L{L+l) 



(16) 



The position of the barrier in Eq. (|T6|) has to be in the 
region where Eq. (fT5|) is satisfied. This happens for suit- 
ably large dipole moments, d > d a where 



d a 



(h 2 L(L + l)f C 6 (4^ £o ) 4 
27^ 3 a 4 



1/8 



(17) 



The subscript a stands for the attractive interaction. For 
two indistinguishable fcrmionic polar 40 K 87 Rb molecules, 
and for L = 1 and Ml = 0, a(l, 1;0) = 4/5 and we get 
d a = 0.103 D. The threshold laws for quenching colli- 
sions in an electric field are the same as in the zero-field 
limit. Consequently, the quenching cross sections and 
rate coefficients behaves as in Eq. (fTT|) except that Vb is 
given now by Eq. (|16p and varies with d. We can insert 
Eq. (|16p in Eq. (fTTj) to get the corresponding analytical 
expressions. For a partial wave L > 0, the quenching rate 
scales as d 4 ^ L+ 2 ) . For two indistinguishable fermionic po- 
lar molecules at ultracold temperatures when L = 1 and 
M L = 0, we get 



K qu 

1x L = 1,Ml=0 



3tt /6V 6 



5 6 h 14 
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47T£r 



k B TxA. (18) 



Thus the L = I, Ml = quenching rate increases as 
d 6 . This is a more rapid dependence on dipole moment 
than for purely long-range dipolar relaxation in dipolar 
gases 



For negative C3 coefficients in Eq. 



-C 3 /i? d 



is re- 



c 3 




h 2 L(L + l) 


R 3 


> 


2/iR 2 



pulsivc in Eq. (JSJ) . For L = 1 partial waves for example, 
this occurs when Ml = ±1, which favors a repulsive ori- 
entation of dipoles. We consider 



(19) 



in Eq. ([8]). The long-range potential again experiences a 
barrier, but now it is generated by the balance between 
the repulsive dipole potential at large i?, and the attrac- 
tive van der Waals potential at somewhat smaller R. In 
this case, the position and height of this barrier are given 
by 



R b 



2C 6 



\C, 



1/3 



4C 6 



<xd 4 



(20) 
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For this approximation to hold, the position of the barrier 
in Eq. (|20| has to be in the region where Eq. (fT9|) is 
satisfied. This requires that d > d r where 



(h 2 L{L + l)f C 6 (4^ ) 4 
4/i 3 a 4 



1/8 



(21) 



The subscript r stands for the repulsive interaction. For 
two indistinguishable fcrmionic polar 40 K 8T Rb molecules, 
and for L = 1 and M L = 1 or M L = -1, a(l, 1;±1) = 
2/5 and we get d r — 0.186 D. We can replace Eq. (|2U|) in 
Eq. (fTTj) to get the corresponding analytical expressions. 
For a partial wave L > 0, the quenching processes scale 
as d~ 4( - i+ 2). For two indistinguishable fermionic polar 
molecules at ultracold temperatures when L = 1 and 
\M L \ = 1, we get 



K qu , 



3tt / 50 h? C R 



3/2 



/' 



4tt£ 



-3 



k B T x A. (22) 



The L = 1, | Mi | = 1 quenching rate decreases as <i~ 6 as 
the electric field grows. 

These analytical expressions use the diabatic barri- 
ers. If we consider that at large d, the total rate is 
mostly given by the Ml = contribution (wc neglect 
the |Ml| = 1 contributions at large d), one can have an 
analytical expression using the adiabatic barrier. If we 
take into account the couplings between L = 1 , Ml = 
and L = 3, Ml = 0, we can diagonalize analytically the 
2x2 matrix in Eq. ([8]). It can be shown that for each 
dipole moment d, the coupling with L = 3, Ml = lower 
the diabatic barrier of L = 1, Ml = by a factor of 0.76 
at the position of the barrier, to give rise to the adia- 
batic barrier. Inserting this correction of the barrier in 
Eq. ([Til , this yields a correction of 0.76 -3 / 2 sa 1.51 for 
:i t M L =o- The difference between diabatic and adia- 
batic calculations can be already seen in Fig. 2] and Fig. [5] 
for the numerical barriers at large dipole moment. 

In Fig. [6] the black curve corresponds to the total 
quenching rate coefficient as a function of d for two in- 
distinguishable fcrmionic polar 40 K 87 Rb molecules in the 
same quantum state for L = 1 and at a temperature of 
T = 350 nK. The analytical expressions I, II, III, IV, 
V (green thin lines) correspond respectively to Eq. ([IB"]) , 
2 x Eq. ([13]), Eq. ([18]), 2 x Eq. ([52]), 1.51 x Eq. (flgf . 
The curves III and IV are for the diabatic barriers, while 
curve V is to account for the adiabatic barrier. The red 
dashed line (I+III) represents the L = 1, Ml = partial 
wave contribution for the diabatic barriers while the blue 
dashed line (II+IV) represents the sum of L = 1, Ml = 1 
and L = 1,Ml = — 1 partial wave contributions for the 
diabatic barriers. The analytical sum I+II+III+IV is 
represented as a black dashed line. To account for the 
adiabatic barriers we assume that the correction for the 
total rate comes only from the L = 1,Ml = partial 
wave, and we replace I+III by I+V (red solid line). The 



analytical sum I+II+V+IV is represented as a black solid 
line. Neglecting the d~ e contribution at larger d, the an- 
alytical p— wave quenching rate (taking into account the 
adiabatic barriers) is given by the simple expression 



3 17 p 3 C^ 1/4 



h 10 



1MP2 



29311^6 X 1/2 

5 6 ft 14 



d 6 



(4^o) s 



k B T x A. (23) 



Pi (jp2) is the quenching probability reached at the height 
of the barrier in the QT model for the zero (non-zero) 
electric field regime. The QT model assumes that p\ = 
P2 = 1 but become fitting parameters {p\,P2 < 1) when 
compared with full quantum calculations or experimental 
data. The limiting value d a = 0.103 D (d r = 0.186 D), 
where the d® (d~ 6 ) behavior begins, has also been indi- 
cated with an arrow. It turns out that the total rates 
for L = 1 calculated analytically (for both the use of 
diabatic and adiabatic barriers) are very similar to the 
numerical ones of Fig. [5] (10 % difference at most, around 
d a )- However, the sub-components L = 1,Ml = and 
L = 1, 1 Ml I = 1 have different behaviors. For example 
the numerical L = l,Mj, = (L = 1, |Ml| = 1) com- 
ponent starts to increase (decrease) at earlier dipole mo- 
ment (typically at 0.02 D) than their analytical analogs 
(typically after 0.06 D). The use of the simple analyt- 
ical expressions (using the diabatic or adiabatic barri- 
ers) can be useful to estimate the total rate coefficients, 
while the numerical ones are prefered to estimate the 
L = 1, M L = and L = 1, \M L \ = 1 individual rates. 



III. PROSPECTS FOR COLLISIONS IN TWO 
DIMENSIONS 

In three dimensional collisions, the quenching loss is 
largely due to incident partial waves with angular mo- 
mentum projection Ml = 0, emphasizing head-to-tail 
orientations of pairs of dipoles. These are the kind of 
collisions that arc largely suppressed in traps with a pan- 
cake geometry, with the dipole polarization axis orthogo- 
nal to the plane of the pancake [37[ . If these collisions can 
be removed, then it is likely that increasing the electric 
field will suppress quenching collisions, making evapora- 
tive cooling possible. If wc assume an ideal pancake trap 
that confines the particles to move strictly on a plane, one 
can apply the present model to estimate the behavior of 
the quenching processes. We assume that the molecules 
are polarized along the electric field axis, perpendicular 
to the two dimensional plane. In this case, the long range 
potential is given by 



V(p) = 



h 2 (\M\ 2 - 1/4) 
2~V 



47T£o p 3 



(24) 



where p stands for the distance between 2 particles in 
a two dimensional plane, M stands for the angular mo- 
mentum projection on the electric field axis. The last 
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0.2 0.3 
d(D) 



trifugal term in Eq. ([8]) leads to the threshold laws in 
Eq. ([10|). the replacement L(L + 1) |M| 2 - 1/4 (that 
is L ^ |M| - 1/2) in Eq. (TO} leads to 



|T A 7| 2 



(25) 



where V;, denotes the height of the centrifugal barrier in 
two dimensions. This result requires that the centrifugal 
potential is repulsive, i.e., that \M\ > 0. For M = 
the threshold law exhibits instead a logarithmic diver- 
gence [3l| . 

In two dimensions, quenching cross sections and rate 
coefficients have respectivelly units of length and length 
squared per unit of time, and are given by (39j 



FIG. 7: (Color online) Barrier heights as a function of the 
induced dipole moment d for the partial waves \M\ = 1 in two 
dimensions. The green thin curves represent the analytical 
Eq. (HU) (constant) and Eq. ([29} (d 4 ). The dashed black curve 
is the sum of them. The solid black curve is the height of the 
barrier in Eq. ([24)) . 



qu 
T M = 



ra 2 x a 



K qu 

M 



- \Tft\ 2 x A. 



(26) 



Within this model, it follows that the quenching cross 
section and rate coefficient for |Af| > are given by 



KRb( Z,v=0,N=0) + KRb( T,v=0,N=0) ; |M|=1 ; T=350 nK 



— analytical limits 

- - total analytical 

— total numerical 




VII 



FIG. 8: (Color online) Two dimensional quenching rate coef- 
ficient (black curves) of two indistinguishable fermionic polar 
40 K 87 Rb molecules as a function of the induced electric dipole 
moment for the M — 1 and M = — 1 components at a tem- 
perature of T = 350 nK. The dashed lines represent the rate 
using analytical expressions while the solid line represents the 
rate using the numerical expression (see text for detail). The 
individual analytical curves have been indicated in the right 
hand side by roman numbers. 



M 



K qu 

M 



2/i V b 
h 



\M\ 



E \M\-l/2 x A 



< k!. M| > xA. 



mv 6 



M\ 



(27) 



The energy dependence is in agreement with the one 
found in Ref. g^j. In Eq. ([53), 



my 



54 n 3 C 6 

for the zero-electric field regime and 
1 / £ 



1/2 



4C 6 \47T£ 



(28) 



(29) 



for the non-zero electric field regime. The height of these 
barriers has been reported in Fig. [7] (green thin lines). 
These results imply that for \M\ = 1 the quenching pro- 
cesses within this model will be independent of the dipole 
moment in the zero electric field regime, where 



K qu 



2VC 6 
ft 4 



1/2 



k B T x A 



(30) 



and will scale as d 4 in the non-zero electric field regime, 
where 



term comes from the repulsive dipolc-dipole interaction 
when the dipoles are pointing along the electric field and 
approach each other side by side. The height of this bar- 
rier has been plotted as a function of d in Fig. [7] (black 
solid line). At ultralow energy and large molecular sep- 
aration, the Bethe-Wigner laws for quenching processes 
depend only on the long-range repulsive centrifugal term 
1/R 2 dl], |38|. The repulsive centrifugal terms are dif- 
ferent in Eq. © and Eq. ([M|) . As the repulsive cen- 



K qu 



AhCe 



d 2 



47T£ 



k B T x A. 



(31) 



We use the fact that < E c >— ksT in two dimensions. 
The non-zero electric field regime is reached when d > 
dir> where 



l 2D 



(ft 2 (|M| 2 - 1/4)) 3 C 6 (4^q) 4 \ 1/8 

4/*3 J 



(32) 
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For two indistinguishable fermionic polar 40 K 87 Rb 
molecules, and for \M\ = 1, we get diu = 0.081 D. 

The behavior of the quenching rate (black lines) is 
shown in Fig. [5] for two indistinguishable fermionic polar 
40 K 87 Rb molecules as a function of the induced electric 
dipolc moment for M = 1 and M = — 1 components at a 
temperature of T = 350 nK. The dashed line represents 
the analytical rate which is the sum of the analytical 
expression VI corresponding to 2x Eq. (|30|) and analyt- 
ical expression VII corresponding to 2x Eq. (f3"Tj) . The 
solid line represents the rate using the general expression 
Eq. (|27[) and the numerical height of the barrier calcu- 
lated in Eq. ([23]). The limiting value d 2 D = 0.081 D, 
where the d~ A behavior for the quenching rate begins, 
has also been indicated with an arrow. The difference 
between the numerical calculation and the analytical ex- 
pression reflects the difference in the calculation of the 
height of the barrier, already seen in Fig. [71 The nu- 
merical calculation is more exact, while the other is an- 
alytical. However, at large d, the numerical rate tends 
to the analytical d~ 4 behavior. The quenching rate de- 
creases rapidly as the dipolc moment increases and this 
may be promising for efficient evaporative cooling of po- 
lar molecules since the elastic rate is expected to grow 
with increasing dipole moment (4lj . 

IV. CONCLUSION 

We have proposed a simple model which combines 
quantum threshold laws and a classical capture model to 



determine analytical expressions of the chemical quench- 
ing cross section and rate coefficient function of 
the collision energy or the temperature. We also pro- 
vide an estimate as a function of the induced electric 
dipole moment d in the presence of an electric field. We 
found that the quenching rates of two ultracold indis- 
tinguishable fermionic polar molecules grows as the sixth 
power of d. For weaker electric field, quenching processes 
are independent of the induced electric dipole moment. 
Prospects for two dimensional collisions have been dis- 
cussed using this model and we predict that the quench- 
ing rate will decrease as the inverse of the fourth power 
of d. This fact may be useful for efficient evaporative 
cooling of polar molecules. This model provides a gen- 
eral and comprehensive picture of ultracold collisions in 
electric fields. Preliminary data suggest that this model 
gives good agreement with experimental chemical rates 
for three dimensional collisions in an electric field 42] . 
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